Stochastic ionization through noble tori : Renormahzation results 
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We find that chaos in the stochastic ionization problem develops through the break-up of a 
sequence of noble tori. In addition to being very accurate, our method of choice, the renormahzation 
map, is ideally suited for analyzing properties at criticality. Our computations of chaos thresholds 
agree closely with the widely used empirical Chirikov criterion. 
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■ The multiphoton ionization of hydrogen in a strong microwave field revolutionized traditional notions about 
the physics of highly excited atoms. Its interpretation remained a puzzle until its stochastic, difFusional nature was 
^ uncovered through the then-new theory of chaos j|] thereby making it the most important testing ground for quantal 
d ' manifestations of classical chaos |^ . The onset of chaos in that problem has been the focus of intensive research and is 
• generally believed to be induced by the break-up of invariant tori Q , which act as barriers in phase space preventing 
[ the diffusion of trajectories for Hamiltonian systems with two degrees of freedom. 
0^ . Over the last two decades, different methods to estimate break-up thresholds of invariant tori have been developed 
[ and applied to various physical models with two effective degrees of freedom: Chirikov's resonance overlap crite- 
I— , i rion [|j , its empirical "2 /3-rule" version , and Escande-Doveil renormahzation ^ . These methods provide merely 
^**^ ' approximate values (as compared to direct numerical integration) even though Escande-Doveil renormahzation is 
known to give very accurate values in some concrete examples. Systematic methods are now available to obtain very 
accurate values, such as, e.g., Greene's residue criterion ||^, Laskar's frequency map analysis |p|,p^, or renormahzation 
analysis [pd|-p^. There is numerical evidence that these independent and systematic methods give the same values 
for the chaos thresholds ||l^,|lj,|l^ . However, renormahzation has an additional advantage : By focusing on specific 
tori, it leads to very accurate thresholds, and allows one to analyze the properties at criticality. Indeed, the renormal- 
ization map can be likened to a phase space microscope by which the system can be studied with larger and larger 
magnification. 

In this article, we find the transition to chaos in the hydrogen atom driven by circularly polarized microwaves 
\^ using renormahzation. This problem has emerged as paradigm for a number of issues in multi-dimensional nonlinear 
dynamics (see Rcf . and references therein) . Our results reveal that the onset of chaos in the microwave problem 
. develops through a sequence of "noble" tori, i.e. tori with frequency equivalent to the golden mean 7 = (\/5 — l)/2. 
' Furthermore, the chaos thresholds obtained by the renormahzation are in very good agreement with the 2/3- rule 
criterion ||l7|-p^, and better than the ones obtained by the Escande-Doveil renormahzation In this model, 

Chirikov's criterion can be used to determine very accurate critical thresholds of ionization (accurate to 5%) for most 
values of the parameter (eccentricity of the initial orbit) . 

Renormalization method- The renormahzation method is based on the construction of successive canonical trans- 
formations 1 11 l^,|ljl. In its most recent version, it acts on the following family of Hamiltonians with two degrees of 



>1 11 II v 
. freedom written in actions A = {Ai, A2) and angles (p = (</?i, (^2) 



■ H{A,ip) =uj ■ A + V{n- A,ip) (1) 

=a;-A+ ^ 14,.(f^-A)'=e^^-'^, 

fceN 

where a; = (lu,—1) is the frequency of the invariant torus, and = (1,0;) is some other vector. Two steps are 
involved [ po[ : an elimination of the non-resonant modes v of the Hamiltonian (the modes which do not involve small 
denominator problems), and a rescaling of phase space (shift of the resonances, rescaling in time and in the actions). 
This transformation reduces to a map TZ of Fourier coefficients {uj', a'; {V^' ^}) = TZ{uj, a; {Vk^i,})- 

The main conjecture of the renormalization approach is that if the torus exists for a given Hamiltonian H, the 
iterates TZ"H of the renormalization map acting on H converge to some integrable Hamiltonian Hq. This conjecture 
is supported by analytical results in the perturbative regime po| , ^ , and by numerical results [p"5|jl4[ . For a one- 
parameter family of Hamiltonians {Hp}-, the critical amplitude of the perturbation Fc{uj) is determined by the 
following conditions : 
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TZ'^Hf Ho{A) =LJ-A+I{n-Af for F < Fc{uj), (2) 
7^"i^F -> oo for F > Fc{lu). (3) 

n — ^C30 

Model- We consider a hydrogen atom interacting with a strong microwave field of amphtude F and frequency Q, 
circularly polarized in the orbital plane. In action-angle variables, the classical Hamiltonian is reduced to the following 
Hamiltonian with two degrees of freedom [|l8| : 

H{J,L,e,iP)^-—-QL + F J2 K(J,i)cos(n0 + ^), (4) 

n— — cxj 

where 

Vo{J,L)^-^J^, 



VniJ,L) = - 

n 



Jnine) + — — —Jnine) 



r, for n 7^0, 



where Jn is the nth Bessel function of the first kind and is its derivative. The angles 9 and V' are conjugate to 
the principal action J and to the angular momentum L respectively. The eccentricity e of the initial orbit is given 
by e (1 - LVj2)1/2. Hamiltonian (|) can be rescaled in order to eliminate the dependence on the frequency of the 
microwave field. We rescale time by a factor J7 [we divide Hamiltonian (Q) by Vt]. We rescale the actions J and L by 
a factor A = Q}/'^, i.e. we replace H{J, L, 9, ip) by XH{J/\, L/X, 9, ip). We notice that this rescaling does not modify 
e{J,L). The resulting Hamiltonian becomes : 

H=-—-L + F' VM,L) cosine + ^P), 

n— — oo 

where F' = Ffi""*/"^ is the rescaled amplitude of the field. In what follows we assume that 1. 

For the unperturbed Hamiltonian i.e. with _F = 0, we consider a motion with Kepler frequency lo G [—1, 1] (the 
high scaled frequency regime). In phase space, this trajectory evolves on a two-dimensional torus. If oj is irrational, 
the trajectory fills the torus densely. This invariant torus is located at J = w"^/'^ and L = —Eq — ^uP'l'^ ^ where i?o is 
the total energy of the system (in the rotating frame). For convenience, we shift the action J such that the torus is 
located at J = 0, and we expand the resulting Hamiltonian in Taylor series in the action J. Furthermore, we rescale 
the actions J and L by a factor A = — 3w^/'^, making Hamiltonian (^) : 

iJ;^=c.J-L + ^-|i^J^- + F(-V4/3j2-H2c.-V3j_3..2/3\ ^ V;(e)cos(n0 + ^), (5) 

where V^(e) — y„(l, Vl — e^). We note that the rescaling coefficient A has been chosen such that the quadratic part 
of the Hamiltonian is J^/2. The resulting Hamiltonian conforms to (0) with a = and with coordinates A = (J, F) 
and y = (9, -0)- 

In what follows, we consider e (the eccentricity of the initial orbit) as a parameter of the system. For e = 1, the 
resulting model is a one-dimensional hydrogen atom in a linearly polarized microwave field By varying e, we 
obtain a wide variety of models where we can compare renormalization and empirical rules. 

Breakup thresholds of invariant tori- For a given eccentricity e, we compute the critical function Fc{uj]e) using 
the renormalization method [Eqs. (§)-(^)]. Figure 1 shows a typical critical function uj i— > Fc(w;e) for e = 0.45 and 
for co-rotating orbits (a; > 0). The critical function vanishes at all rational values of the frequency (since all tori with 
rational frequency are broken as soon as the field is turned on), and it is discontinuous on a dense set. This figure shows 
that there is no invariant torus between the resonances 1:1 and 4:1 for Fc > max^g[i/4 ij Fc{uj;e = 0.45) ~ 0.0069. 
The most stable region is located around oj « 0.723. 

Figure 2 shows the critical thresholds 

between two primary resonances maX(^g[i/2,i] Fc{uj,e) and maXj^g[i/3 1/2] ^c('^,e) as a function of the parameter e 



for CO- rotating (w > 0) and for counter- rotating (w < 0) orbits. This figure is analogous to Fig. 3 of Ref. |19|. In a 



broad range of values of e (e e [0.2, 1]) the heuristic 2/3-rule criterion gives very accurate results (accurate to 5%) by 
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comparison with our renormalization results for co-rotating orbits. The values are even better than the ones computed 
by Escande-Doveil renormalization JTof . However, for very low eccentricity, there is a large gap between both results 
which makes the 2/3-rule inapplicable in this range of parameter e. For this case, Escande-Doveil renormalization 
is a better criterion for the determination of the thresholds. As e tends to zero, the discrepancy is even bigger. For 
instance, for the overlap between 1 : 1 and 2:1, the 2/3-rule criterion predicts a finite value of the threshold at e = 
whereas Hamiltonian (^) is integrable in that case and the critical function is expected to go to infinity at e = 0. 

For counter- rotating orbits, the 2/3-rule criterion overestimates the critical couplings even though it gives fairly 
accurate results (less than 10% for e G [0.8,1]). The resulting critical curve is below what has been obtained in 
Ref. using the Escande-Doveil renormalization with a discrepancy around 30%. 

Our results confirm that the orbits with medium eccentricity can diffuse more easily between the first primary 
resonances (n = 1,2,3) However, in order to ionize, the orbits must diffuse throughout phase space or at least 
between a large number of primary resonances (in experiments and numerical simulations p^ , « 40). Since 
the 2/3-rule gives very accurate results in a broad range of parameter e, we compute the critical thresholds between 
resonances m:l and m-|-l:l for m = 1, . . . ,n1 with = 40, for medium and large eccentricities. In Fig. 3 we have 
plotted the different critical curves e i-^ maX(jg[i/(„j+i) i/,„] Fc{lu; e) for m — 1, . . . , for co-rotating and counter- 
rotating orbits. It appears that there is a very broad region of the parameter e where an orbit cannot diffuse from the 
resonance 1:1 to n^:l. In fact, only in the region e S [0.9, 1] (with F > 0.015 for co-rotating orbits, or with F > 0.05 
for counter-rotating ones), the orbits can ionize. In other terms, only the high-eccentricity orbits can ionize in this 
classical model (||). This observation reinforces the importance of core collisions for the ionization process [p^ . 

A more commonly used function is the scaled function Fo{uj) = Fc{uj)uj^'^^'^ . For the linearly polarized case (e = 1), 
experiments |23) and numerical computations j2j] show that for Fq = FcTi^ > 0.02 (when n is the principal quantum 
number) there is ionization for fio = 1/w (the scaled frequency) close to 1. Here we find agreement with this value: 
for _Fo > 0.022, there is no invariant curve in phase space and diffusion can occur. For the circularly polarized case 
(e 7^ 1), the conclusions are less clear. However since there is a broad stable region e e [0, 0.8], we expect the circularly 
polarized driven atoms to ionize less easily than the linearly driven ones. Since the orbits of high eccentricities are 
the ones which ionize, we expect the behavior of ionization curves for circularly polarized microwaves to be similar to 
the ones for linearly polarized ones. This is consistent with experiments (Fig. 1 of Ref. j2^ ) . 

Progress to chaos through noble tori- The renormalization map allows us to determine accurately the frequency of 
the last invariant torus to break-up. The conventional wisdom is that the last invariant torus surviving with increasing 
the amplitude of the perturbation is the one with frequency equal to the golden mean 7 = {V5 — l)/2. The example 
we study shows that this belief is mistaken. Figure 4 shows the value of the frequency of the last invariant torus 
between resonances 1:1 and 2:1, and between 2:1 and 3:1 (for co-rotating and counter-rotating orbits), as a function 
of the parameter e. We have identified the frequency of these tori (by accurate computation of critical couplings in 
the neighborhood of these frequencies). In the range of the parameter e and for these regions of phase space, each 
last invariant torus is a noble one in the sense that its frequency lu is equivalent to the golden mean 7, i.e. there exist 
integers (a, 6, c, d) g such that ad — be = ±1 and u — (aj + b)/{ej + d) (or equivalently the tail of the continued 
fraction expansion of w is a sequence of 1). For example, for e = 0.45, the last invariant torus is expected to be 
(7 -I- l)/(27 + 2) « 0.7236 (see Fig. 1). For the region of high eccentricities (which ionize more easily), the expected 
frequency for the last invariant torus is = (27 + l)/(37 + 2) « 0.7835. The observation that the last invariant 
torus is noble has also been checked for the counter-rotating case. Whether or not this observation holds for a generic 
Hamiltonian remains an open question [^5| . 

Conelusion- Using the renormalization method, we find that the empirical 2/3-rule Chirikov criterion is surpriz- 
ingly accurate for the onset of chaos in the stochastic ionization problem. The model studied in this Letter, the 
hydrogen atom driven by strong microwaves, shows how renormalization and empirical rules can be used together in 
order to obtain very accurate information on the stability of the system. 
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FIG. 1. Critical function Fc{uj) for lu G [0, 1] for Hamiltonian with e = 0.45. 
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FIG. 3. Critical threshold Fc between various primary resonances m:l and m+l:l for m = 1, . . . , 40 from left to right, 
obtained by the 2/3-rule criterion for co-rotating orbits (continuous curves) and counter-rotating orbits (dashed curves). 
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FIG. 4. Value uJiast of the frequency of the last invariant torus for Hamiltonian as a function of the parameter 
between resonances 1:1 and 2:1, and (6) between resonances 2:1 and 3:1. 
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